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Abstract 

We present a first principle study of new class of high-T c half-heusler ferromagnets NiCrZ (Z = 
Si, P, Ge, As, Te). The structure and magnetic properties are investigated through the calculation 
of the electronic structure, equilibrium lattice constant, magnetic exchange interaction Jij and 
Curie temperature T c . The role of sp-elements and the influence of lattice expansion/compression 
are also studied. In alloys having 20 valence electrons, a pseudo-gap of the majority band can be 
formed at Fermi level. Otherwise, the half-metallicity and ferromagnetism at temperatures much 
higher than room temperature are found to be stable in a wide range of lattice expansion. Based 
on these results, NiCrZ can be expected to be promising materials for spintronics. 
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To exploit the great potential of spintronics, seeking for new magnetic materials and 
half metallic ferromagnets is one of the most important issues in materials physics. Half 
metallic ferromagnets, whose the electronic structure behaves like metals w.r.t the electrons 
of one spin direction and like semiconductors w.r.t. the spins in the opposite direction, 
have the extreme spin polarization at Fermi level Ep and can achieve the ferromagnetism 
at temperatures higher than room temperature. 

Since the pioneer work on the half-metallicity in half-heusler alloy NiMnSb by de Groot 
and collaborators |l| in 1983, heusler and half-heusler alloys have been attracted great 
attention in both theoretical and experimental study. Especially, the family of half-heusler 
alloy XYZ (X and Y are transition metals, and Z is an sp— valence element), has been one 
of the most interested compounds , [§]] due to the possibility of half-metallic 

: erromagnetism being achieved. Besides the extreme spin polarization has been verified 
4], [5I, 0], the ferromagnetism in the half- Heusler NiMnZ at temperatures ranging between 
500K and 730K has been reported for Z= Pd, Pt and Sb [2] by experiments. Theoretically, 
the first principle study has predicted the half-metallic character in their electronic structure 
of many half-heusler alloys such as FeMnSb, CoMnSb and FeCrSb, etc. js, 7]. Recently, the 
half-metallic ferromagnetism with very high Curie temperature (T c ) in a new class of half- 
heusler alloys NiMnZ (particularly, T c m 1000K for NiMnSi) has been predicted Also, 
the half-metallic ferromagnetism with T c higher than room temperature can be expected 
even if Cr replace Mn in NiMnZ. 

In this Letter, based on the first principle study of the structural and magnetic proper- 
ties of NiCrZ (Z = Si, P, Ge, As and Te) we propose a new class of high-T c half-metallic 
ferromagnets. For this aim, we have carried out the calculation as follows. First, we perform 
the total energy calculation to evaluate the equilibrium lattice constant (ELC) by means of 
the generalized gradient approximation (GGA) within the framework of ultrasoft pseudo- 
potential (UPP) (thanks to STATE-Senri) and both of the muffin-tin (MTA) and atomic 
sphere approximations (ASA) within KKR-LSDA (thanks to AKaiKKR). Using ELC, we 
next calculate the density of states (DOS), magnetic exchange interaction between mag- 
netic sites. Finally, using as the input data we calculate T c by employing three statistical 
approaches: mean field approximation (MFA), random phase approximation (RPA) and 
Monte Carlo simulation (MC). 

Half-heusler alloys NiCrZ have C1& structure. The study of structural property 0, [§} 
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TABLE I: ELC (au) by means of GGA within UPP (a UP ), MTA (o M t) and ASA (a A s) within 
KKR-LSDA. 



Alloy 


a\jp 


a AS a MT 


NiCrSr 


10.2714 10.3243 10.4450 


NiCrP 


10.2622 


10.3823 10.5390 


NiCrGe' 


10.4521 


10.5523 10.6652 


NiCrAs' 


10.5631 


10.6310 10.7966 


NiCrTe' 




11.5116 11.7980 



shows the most stable configuration of NiYZ should be in a phase (except Z being a lighter 
sp-element such as C and N) with Ni being located in the octahedral coordinated pocket 
(7,7,4), ^ a ^ an sp-element Z at the body center (|, 1) anc ^ a vacan t site at 

(|, |, |). Ni atom possesses 10 valence electrons and Cr has 6 electrons. Therefore, the 
numbers N t of valence electrons per formula unit are 20, 21 and 22 corresponding to sp- 
elements of group IV, V and VI, respectively. Regarding to "rule of 18" of the half-heusler 
alloys M = N t — 18, the total magnetic moment M of NiCrZ should be 2//_b, 3/jb and 4//g 
corresponding to Z of the groups IV, V and VI, respectively. 

Table [I] shows ELC calculated by means of three approximations: GGA within framework 
of UPP, MTA and ASA within LSDA. As can be seen from Tab. QJ MFA gives the lattice 
constant aMT larger than ASA and UPP, whereas UPP gives the smallest one. However, the 
deviation between aup, o>as and a MT is considerably small (a^s — o-up ~ 1% f° r NiCrSi). In 
general, the lattice constant is larger if the atom radius is larger. However, UPP gives aup 
of NiCrP smaller than NiCrSi. Referring to the case of NiMnZ js], we can also expect that 
the real lattice constant of NiCrZ might be arranged in the range from aMT to aup- 

The left panel of Fig. [1] illustrates DOS of the typical alloys standing in the compounds 
which possess the different numbers of valence electrons. The DOS of three alloys possess- 
ing 20 valence electrons (NiCrSi), 21 valence electrons (NiCrAs) and 22 valence electrons 
(NiCrTe) are plotted. The change of the electronic structure w.r.t the lattice constant is 
shown in the right panel for NiCrAs as a representative case. In addition, it should be noted 
that NiCrZ alloys with Z = Sn, Sb, Se and S are metals at ELC. Furthermore, NiCrZ with 
Z being an element of group IV exhibits a pseudo-gap of the majority band at Ep in UPP 
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FIG. 1: DOS of NiCrSi, NiCrAs and NiCrTe at a = clas are given in left panel. Volume dependence 
of DOS within LSDA and GGA is shown for NiCrAs in right panel . 

calculation, and its electronic structure shows the low spin behavior at ELC but the high 
spin picture at sufficiently larger a (about 0.5% larger than ELC). 

As seen from the left panel of Fig. [TJ while almost states of 3d states of Ni concentrate 
at lower energies (in valence band) and are occupied, the 3d states of Cr are distributed at 
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both of the lower and higher energies and strongly polarized. The 3d— majority spin states 
of Cr distribute at the lower energies and are almost occupied. These states hybridize with 
3d states of Ni (and 3p of Z elements) to form a band at Ep. The minority spin states of 
Cr shift to the higher energies and are unoccupied, leading to the gap (that Ep falls into) 
being formed in the minority band and the half-metallic character. For alloys possessing 
21 valence electrons, Mn-3d states in the majority gap split into two bands e g and t a , and 
Ep lies at the vary narrow valley which is formed between these bands. Moreover, while 
the half-metallic character is obvious in DOS calculated by LSDA, GGA combined in UPP 
gives a very narrow pseudo-gap at Ef in the majority band of NiCrGe and NiCrAs, and 
the half-metallic behavior might be disappeared. The calculated total magnetic moment M 
of NiCrSi is exactly in consistence with "rule of 18". Similarly, we have also obtained 
the half-metallic behavior of the remaining alloys NiCrZ (Z = P, Ge, As, Te) and the total 
magnetic moment of 2fj, B for Z of group IV and 3\xb for Z of group V at a = aup, a as an d 
dMT (except NiCrGe and NiCrAs at a = clmt- See also Tab. HI!), an d 4/ie for Te at a = clas 
and (Xmt ■ 

To study the influence of the lattice expansion (or compression), we draw the DOS of 
NiCrAs at three different ELC in the right panel of Fig. [IJ LSDA (black line) and GGA 
(blue line) calculations are given. Since the minority states of Cr-3(i electrons in GGA result 
is shifted toward the higher energies while the majority band is mostly kept unchanged, the 
minority gap in GGA is wider than LSDA. For NiCrAs, the half-metallicity is destroyed at 
omt- As the lattice constant a decreases, Ep shifts toward the conduction band, the half- 
metallicity becomes more stabilized; however, the half-metallicity can be destroyed if a is 
smaller than the threshold value at which Ep falls into the conduction band. Except NiCrAs 
and NiCrGe whose the half-metallicity might be destroyed at a > cimt-, our calculation shows 
that the proposed alloys preserve the half-metallicity in a wider range of the expansion 
(compression) of the lattice cells. 

To investigate the ferromagnetism in NiCrZ alloys, we calculate the magnetic exchange 
interaction J{j at three values of ELC. The exchange interaction Jy between two impurities at 
ith and jth sites, which are embedded in the ferromagnetic medium, is efficiently calculated 
by utilizing the magnetic force theorem. The frozen potential approximation [9|] is employed 
and Liechtenstein formula [lOjj is used for evaluation of Jy. J^- of NiCrZ at a = a as is 
illustrated in FigJ21 For alloys with N t = 20, the magnetic exchange interaction is 
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FIG. 2: Cr-Cr effective exchange coupling constant Jjj vs. distance d in units of lattice constants 
of NiCrZ. Jij is calculated at lattice constant a = a as- 



ferromagnetic at the 1st, 3rd and 5th nearest neighbors, but anti-ferromagnetic for the 2nd 
and 4th nearest neighbor pairs. For N t = 21, is ferromagnetic for the 1st, 2nd 5th and 
m nearest neighbor pan,. Similar to NiMnZ Q, J„ becomes ferromagnetic at distances np 
to the 5th nearest neighbors for N t = 22. Thus, Jo? with j — 1 and j = 2 becomes larger 
with increasing N t . For NiCrTe, due to the larger distance between atoms, Joi becomes 
remarkably smaller than that of the other half-heusler alloys. 

In order to study the volume expansion effect on the magnetic exchange interaction as 
well as T c , we demonstrate the change of Jy versus the lattice expansion in Fig. El 
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FIG. 3: Jij vs. lattice expansion. 



of the typical half-heusler alloys corresponding to N t = 20, 21 and 22 is calculated for the 
nearest neighbor pairs at lattice constants aup (circles), a as (squares) and clmt (triangles). 
Generally, J02 of all alloys decreases with increasing a. J01 of alloys possessing 21 valence 
electrons increases with a, whereas J01 of alloys with N t = 20 and 22 decreases. Since the 
dominant contributions to T c come mostly from the nearest neighbors at short distances 
(say, at the 1st and 2nd neighbor pairs), T c of alloys with Nt = 20 and 22 actually decreases 
with a. For NiCrP, T c can slightly increases in the range from aup to a as because the J01 
considerably increases, but since Joi's change according to an increase of a at the range of 
larger a is very small whereas J04 becomes more anti-ferromagnetic, T c decreases for a > o,as- 
After calculating J^, the calculation of T c is carried out by three statistical approaches: 
the Monte Carlo simulation T C MC , the mean field approximation (k^ T^ = | J2j Joji where 
ks is Boltzmann constant), and the random phase approximation fill . \\\ ((fc^T^) -1 = 
(^(0) — where J(q) is the Fourier transform of the exchange parameter J«). 



3 

2jv 



Using J^, we evaluate T c . To perform MC simulation, the Metropolis algorithm 13j is 



applied to calculate the thermal average of the ma gne tization M and its powers. Then, 



the cumulant crossing method proposed by Binder 



13| is employed and the fourth order 



cumulant U4 is calculated as a function of temperature for different cell sizes (14 x 14 x 



TABLE II: Predicted T c (K) calculated by MFA, RPA and MC at a = a M T,a A s, aup- 



Alloy 


a 


rpMF 


rpRP rpMC 


NiCrSi 


fl TT T) 


778 


606 


634 


(Nt = 20) 




758 


591 


615 






605 


477 


493 


NiCrGe 


11 TT D 


669 


562 


550 


(A, = 20) 




587 


465 

t: U "J 


47D 


NiCrP 


/Inn 


974 


777 


805 


fiV/ = 21) 




992 


784 


810 




Qj A/fT 


953 


739 


782 


NiCrAs 


«[/P 


1044 


808 


858 


(iV* = 21) 


aA5 


1025 


785 


829 


NiCrTe 


OA5 


665 


588 


605 


(N t = 22) 


aMT 


601 


544 


562 



14, 16 x 16 x 16, and 18 x 18 x 18 conventional fee cells) to find the universal fixed-point at 
T c . Obtained results are shown in Tab. [Til 

It should be emphasized that with a similar method our theoretical T^f c (745K) of 
NiMnSb js| is in very good agreement with the experimental result (T c = 730K) [2]. Similar 
to the case of NiMnZ 8j], MFA overestimates T c even at 100% magnetic atoms, whereas RPA 
underestimates. As usual, the MC simulation gives the T c between values obtained by RPA 
and MFA at the same lattice constant. As expected, T C MC of all alloys is higher than room 
temperature and arranges from 858K to 550K (at ). The lowest T c corresponds to 

NiCrGe with T c = 550K at a = ajjp- The highest T c (858K at ajjp) corresponds to NiCrAs 
and the next are NiCrP (805K), NiCrSi (634K) and NiCrTe (605K). It is noted that, obtained 
T c here is calculated based on the electronic structure within KKR-LSDA. GGA combined 
in UPP shows a pseudo-gap in the majority band of NiCrZ which possesses 20 valence 
electrons; hence, UPP calculation will shows a lower T c for NiCrGe and NiCrSi. Except 
NiCrAs and NiCrGe whose the half-metallicity is considerably sensitive with an increase of 
the lattice constant because of E F closed to the valence band edge, the half-metallicity in 
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the remaining alloys can be preserved when a varies in the range of aup < a < aur- 

In summary, we have investigated the structure and magnetic properties of the half- 
heusler alloys NiCrZ (Z = Si, P, Ge, As, Te). The phase stability is considered. The ELCs 
are predicted by three approximations within the pseudo-potential method and KKR-LSDA 
as well. Using the obtained equilibrium lattice constants, the magnetic exchange interaction 
is calculated and then the evaluation of T c is performed by employing MFA, RPA and Monte 
Carlo simulation. The volume dependence of the electronic structure and magnetic exchange 
interaction as well as T c are also discussed. Following the obtained results, we propose a new 
half-metallic high-T c half-heusler alloys NiCrZ (Z = Si, P, Ge, Te) whose half-metallicity can 
be preserved in a wide range of lattice expansion. 
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